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We present an experimental and theoretical investigation of two-photon direct frequency-comb 
spectroscopy performed through velocity-selective excitation. In particular, we explore the effect of 
repetition rate on the 5Si/2 — > 5D 3 / 2 ,5/2 two-photon transitions excited in a rubidium atomic vapor 
cell. The transitions occur via step-wise excitation through the 5Pi/2,3/2 states by use of the direct 
output of an optical frequency comb. Experiments were performed with two different frequency 
combs, one with a repetition rate of « 925 MHz and one with a repetition rate of w 250 MHz. The 
experimental spectra are compared to each other and to a theoretical model. 

PACS numbers: 42.62.Fi, 32.70.-n 



I. INTRODUCTION 

The use of ultra-short pulses to perform precision spec- 
troscopy has gained increased interest due to the de- 
velopment of phase-stabilized optical frequency combs 
(see Ref. [1] for a review of direct frequency-comb spec- 
troscopy). One of the chief advantages of using the 
output of an optical frequency comb to directly excite 
atomic transitions is the spectral versatility afforded by 
the comb. Stabilized optical frequency combs of the types 
used here contain ~ 10 5 — 10 6 optical frequencies span- 
ning the entire visible spectrum. Thus, the optical fre- 
quency comb acts effectively as a large number of cw 
lasers with frequencies given by 

cj„ = 2w(nf r + f ), (1) 

where n is an integer denoting the mode of the comb, 
f r is the repetition rate of the laser, and fo is the 
carrier-envelope-offset frequency (for a review of fre- 
quency combs see, e.g., Ref. 0). 

Two-photon direct frequency-comb spectroscopy with 
a fully stabilized comb was first applied to the study of 
cold rubidium atoms [3H5]. In these experiments inde- 
pendent control of fo and f r was used to explore the 
effects of resonant enhancement of the two-photon tran- 
sition rate due to resonance with an intermediate state. 
Two-photon direct frequency-comb spectroscopy has also 
been applied to the study of room-temperature atomic 
cesium confined to a vapor cell [5] . The transitions stud- 
ied in that work were excited via two-photon excita- 



tion through a resonant intermediate state using velocity- 
selective excitation. 

Here, we further explore velocity-selective, two-photon 
excitation both theoretically and experimentally. In par- 
ticular, we consider the effect of the repetition rate of 
the comb on the spectra. We identify two different 
regimes: one in which the repetition rate of the comb 
is larger than the Doppler width of the resonance cor- 
responding to the first stage of the transition and one 
in which the repetition rate is smaller than the Doppler 
width of the resonance of the first stage of the transi- 
tion. We present data for the 5Si/ 2 — > 5Pi/ 2 — > 5D 3 / 2 
and 5Si/2 ~~ * 5P 3 / 2 — > 5D 3 / 2 5 / 2 transitions in atomic ru- 
bidium taken with two different optical frequency combs 
with different repetition rates. We compare our results 
to a model calculation and see how the velocity selection 
process differs for the two cases. In addition, we explore 
how the energy of the intermediate state affects the two- 
photon spectra and explain the qualitative differences in 
the spectra that arise from excitation through different 
intermediate states. 



II. THEORETICAL CONSIDERATIONS 

A. Velocity Selective Resonance 

The two-photon transition probability between two 
states I g) and |/) is given by [7] 
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where v is the velocity of the atom, d is the electric dipole 
operator, e-i( 2 ) is the unit vector along the direction of the 
polarization for the first(second) laser beam of intensity 
-^1(2) j Ik and 7/ are the homogeneous line widths of the 
states \k) and |/), respectively, Lo g .^ and ui g: f are the res- 
onant angular frequencies of the transitions \g) — > \k) 
and \g) — > |/}, respectively. The sum over k includes all 
possible intermediate states. We have assumed that the 
laser frequency ui is near the single-photon resonance of 
the \g) — > \k) transitions and have neglected the possibil- 
ity of ll>2 exciting the first stage of the transition. If we 
further assume that there is a single intermediate state 
that is close in energy to o>i we can reduce the sum to a 
single term consisting of the nearly resonant intermedi- 
ate state. The maximum transition probability will occur 
when both the two-photon resonance condition and the 
single-photon resonance condition are met: 



ujg.f =uj\ + ki • v + uj 2 

OJg:k =^1 + kl • V. 
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When an optical frequency comb is used to excite an 
atomic transition, the transition rate must be summed 
over all of the optical modes in the comb [the frequencies 
of which are given by Eq. Q]. Resonance will occur if 
there are any comb modes n\ and n 2 resulting in frequen- 
cies uii and u> 2 that simultaneously satisfy Eqs. ([3| and 
Q. If we consider the situation where the comb light is 
counter-propagated through the atomic sample (so that 
ki = — k 2 ) we can write the resonance conditions as 



-^L =(m + n 2 ) f r + (m - n 2 )— f r + 2/ 
Zir c 



~2tT 



=(»l/r + /o) 1 + 



C 



(5) 
(6) 



where we have defined the x axis to be the direction of 
propagation of the light and taken v x to be the x com- 
ponent of the atomic velocity. For a given n\, n 2 , and 
fo there exists a repetition rate, f r , which will satisfy 
the above equations for some value of v x . If the veloc- 
ity class corresponding to v x that satisfies the resonance 
conditions is populated in the atomic distribution, there 
will be resonant excitation to the state |/) in the atomic 
sample. 

In this work we have investigated this velocity-selective 
excitation with two different frequency combs having dif- 
ferent repetition rates. The first situation we consider is 



the case where the repetition rate of the comb is larger 
than the Doppler width of the resonance corresponding 
to the first stage of the excitation. In this case there will 
exist, at most, one comb mode n\ which is resonant with 
a velocity class present in the atomic sample. As a result, 
there will be at most a single pair of comb modes n\ and 
n 2 that sum up to a total mode number nx = ni + n 2 
that satisfy the resonance conditions given by Eqs. ^ 
and (JsJ) for a velocity class present in the atomic sam- 
ple. As the repetition rate is scanned the velocity class 
excited to the first stage of the transition changes as the 
frequency ui sweeps over the Doppler broadened transi- 
tion. Once the velocity class that is excited to the inter- 
mediate state |fc) also satisfies the resonance condition 
for lu 2 — 2Tr(n 2 f r + fo), the atoms will be excited to the 
final state |/). This results in a sub-Doppler resonant 
peak. 

If the repetition rate is increased over a large enough 
range a different pair of comb modes, adding up to n' T = 
ni +n 2 — 1 will satisfy the resonance requirement and the 
resonant peak will recur. However, because the velocity 
class that satisfies the resonant condition will, in general, 
be different, the relative amplitude of the n' T peak may 
vary greatly as compared to the nx peak; reflecting the 
relative population in the resonant velocity classes. The 
amount the repetition rate must scan before the spectrum 
repeats is 
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If the laser beams are exactly counter-propagating 
throughout the excitation region then the width of the 
resonance in terms of the repetition rate, 8f r , is approxi- 
mately the natural line width of the excited state divided 
by the total mode number: 
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The second case we consider is that in which the repeti- 
tion rate of the frequency comb is less than the Doppler 
width of the resonance corresponding to the first stage 
of the transition. In this case there are multiple comb 
modes resonant with the first stage of the transition for 
different velocity classes. Consequently, there are multi- 
ple pairs of comb modes (ni, n 2 ), having the same total 
mode number tit — n\ + n 2 , that will lead to resonant 
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excitation for different velocity classes. In general, the 
different velocity classes will be resonant at different rep- 
etition rates. The result is a spectrum that has a comb- 
like structure as different pairs of comb modes contribute. 
Figure [T] (a) shows the excitation of the different veloc- 
ity classes to the intermediate and final states for three 
different values of the repetition rate. The resonant repe- 
tition rate, f r .o, is the repetition rate at which the largest 
number of atoms are are excited to the final state. This 
corresponds to excitation of a velocity class that is near 
zero. There are additional resonances when the repeti- 
tion rate is below (above) the resonance value, leading 
to excitation of velocity classes that have negative (posi- 
tive) velocity along the propagation direction of the laser 
beam that excites the first stage of the transition. 
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FIG. 1: Population of the intermediate state and excited state 
as a function of velocity for three different values of the rep- 
etition rate. The bottom trace shows the ground state veloc- 
ity distribution for a rubidium vapor cell at a temperature 
of T = 373 K. As the repetition rate is varied the velocity 
classes excited to the intermediate state move through the 
distribution. When the repetition rate satisfies the resonant 
condition for a given (n\, ri2) and for a specific velocity class 
atoms are excited to the final state. Figure (a) shows the case 
where the resonant repetition rate is different for the different 
mode pairs and velocity classes, leading to multiple peaks in 
the spectrum. The calculation is based on a three-level model 
where the frequencies of the intermediate and final states cor- 
respond to the frequencies of the 5D3/2 state excited through 
the 5Pi/2 state. Figure (b) shows the case where the interme- 
diate state is close to the halfway point of the energy between 
the ground and final states. In this case, the resonant repeti- 
tion rate will be the same for the different velocity classes and 
different mode numbers that have the same tit = ni+n2- The 
frequencies used in the calculation correspond to the frequen- 
cies of the 5D5/2 state excited through the 5P3/2 state. The 
parameters used for the comb light correspond to those used 
in the experiment with the smaller repetition rate frequency 
comb {fr « 250 MHz). 



The case where the repetition rate of the comb is less 
than the Doppler width of the first stage of the tran- 
sition exhibits a significant change if the energy of the 
intermediate state is nearly halfway between the final 
and ground states' energy difference. This is the case 
for the 5Si/2 — > 5P3/2 — > 5D 5 / 2 transition in rubidium, 
where the intermediate state is only f.5 THz from the 
mid-point of the ground-to-final-state energy separation. 
In this case, the different velocity classes satisfy the res- 
onance conditions at nearly the same repetition rates. 
This can be understood by noting that the Doppler shift 
is wavelength dependent. As a result, the apparent spac- 
ing of comb modes in "velocity space" is similar if the 
wavelengths for the first and second stages of the transi- 
tion are similar and are different for the two stages if the 
wavelengths are different. In the case of where the wave- 
lengths are similar the spectrum is significantly simpler 
than the more general case described above and there is 
no comb-like structure in the spectrum. This situation 
is illustrated in Fig. [I] (b), where the energy levels were 
based on the 5Si/2 — > 5P3/2 — > 5D 5 / 2 transition in ru- 
bidium. The resonant peak contains contributions from 
multiple velocity classes that are simultaneously excited. 
Figure [2] shows the calculated fluorescence spectra as the 
repetition rate is scanned for the three-level model both 
for the case where the energy of the intermediate state is 
nearly half that of the two-photon energy (1.5 THz away) 
and for the case where the intermediate state is not close 
to the half of the two-photon energy (7.5 THz away). 
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This dependence of the character of the spectrum on 
the intermediate state energy can be understood by con- 
sidering two pairs of modes (m, n 2 ) and (n[, n' 2 ) = 
(ni + 1, n% — 1) that sum up to the same total mode 
number tit- We suppose that f r is the resonant rep- 
etition rate for a velocity class v x with mode numbers 
(ni, n 2 ) and f' r is the resonant repetition rate for a dif- 
ferent velocity class v' x excited by frequencies with mode 
numbers (n[, n' 2 ). The condition that the two resonant 
peaks overlap is that A/ r = /,', — f r < 2 J» r ■ Solving the 
resonance conditions [Eqs. ^ and ^] gives 
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where we have kept only terms first order in Af r and 
have assumed f r 3> fa for simplicity. In the case of the 



excitation of the 5S! 



/2 



5P 



3/2 



5D 5 /2 transition with 



a comb that has f r ~ 250 MHz we have ri\ k, 1.536 x I0 6 
and n 2 m 1.544 x 10 6 , giving Af r « 0.4 Hz. This is 
smaller than the resonant line width of Sf r ss 0.9, where 
we have used V- ~ 3 MHz for the homogeneous line width 
of the transition. 
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B. Two-Photon Transitions in Rubidium 
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FIG. 2: Calculated fluorescence spectra for a three-level sys- 
tem. The top plot shows the calculated spectrum for the situ- 
ation where the intermediate state is close to half of the two- 
photon transition energy (1.5 THz from the mid point). The 
energies used correspond to those of the 5Si/2 — > 5P3/2 — ► 
5D 5 /2 transition. The different velocity classes are all reso- 
nant near the same repetition rate, resulting in a single peak, 
consisting of unresolved contributions from different veloc- 
ity classes. The bottom plot shows the calculated spectrum 
for the situation where the intermediate state is far from 
the mid point of the two-photon transition energy (7.5 THz 
from the mid point). The energies used correspond to those 



of the 5Sj 
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transition. The param- 



eters used for the comb light correspond to those used in 
the experiment with the small repetition rate frequency comb 
(f r w 250 MHz). 



Figure [3] shows an energy level diagram for the two 
stable isotopes of rubidium ( 85 Rb abundance 72.2% and 
87 Rb abundance 27.8%). The transitions considered here 
are transitions to the 5D 5 / 2 and 5D 3 / 2 states through 
the 5Pi/2 and 5P 3 / 2 intermediate states. The nonzero 
nuclear spin of both isotopes (7 = 5/2 for 85 Rb and I = 
3/2 for 87 Rb) result in hyperfme structure for all of the 
states involved. Thus, there are a number of possible 
transitions leading to a complex and rich fluorescence 
spectrum. Including the hyperfme structure, the general 
two-photon transition rate becomes 
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where Mp, M' F , and M' F are the projections of the to- 
tal angular momenta F, F' , and F" along the axis of 
quantization, ^nh,, is the homogeneous line width of the 



state \nLj), and 0J n LjF:n' V f> is the resonant angu- 
lar frequency of the transition \nLjF) — > \n' L' jr F'). 
The sum over J' runs from 1/2 to 3/2, and the F' runs 
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over the hyperfine levels of the intermediate state. For 
given polarizations of the two beams, the Wigner-Eckart 
theorem can be used to relate the matrix elements to a 
reduced matrix element that is independent of the mag- 
netic sublevels. The reduced matrix elements in the F 
basis can be related to the reduced matrix elements in 
the J basis using standard angular momentum relations 
(see, e.g., [H]). Using these relations, the spectra were 
calculated by inputting known values for the transition 
frequencies [9HI2] and integrating the transition proba- 
bility over the velocity distribution. The number of comb 
modes used for the calculation varied depending on the 
repetition rate being considered. Typically, 15-20 modes 
were used in the calculations that were done for the high 
repetition rate laser (f r w 924 MHz) and 20-25 modes 
were used for the calculations for the lower repetition 
rate laser (/ r « 250 MHz). 

The calculation also accounted for depletion of the 
power for the first stage of the transition as the beam 
propagated through the vapor cell. The two-photon tran- 
sition rate is proportional to the product of the intensities 
of the resonant laser fields. As a result it is advantageous 
to focus the light in order to increase the transition rate. 
This results in a localized region that provides the largest 
contribution to the fluorescence signal. However, if the 
atomic density is large enough, the power resonant with 
the first stage of the transition can be depleted through 
absorption before the focused region of excitation. Since 
the different resonances correspond to excitation of the 
different velocity classes this power depletion varies for 
the different resonance peaks. This effect was taken into 
account in the calculation. For a given resonant veloc- 
ity class, v, the number of atoms resonant with the first 
stage of the transition was calculated using a Voigt profile 
based on the Doppler distribution and number density 
derived from the saturated vapor pressure. The deple- 
tion of the resonant intensity for the first stage of the 
transition was then modeled as 

h = I w e- a °t n{v ^ L , (11) 

where, Jio is the undepleted intensity, cr is the resonant 
scattering cross section for the transition of interest, £ is 
the isotopic abundance, n(v) is the density of atoms in 
the resonant velocity class, and L is an effective scatter- 
ing length. The effective scattering length was adjusted 
to provide good qualitative agreement with the experi- 
mental data. Power depletion of the light resonant with 
the second stage of the transition is unlikely to be sig- 
nificant since the population of the intermediate state is 
much less than the population of the ground state. 



III. EXPERIMENTAL SETUP 

Data were collected by use of two different experimen- 
tal setups. One of them utilized a frequency comb based 
on a TkSapphire laser located at Oberlin College. The 




FIG. 3: Energy level diagram for Rb showing the relevant 
energy levels. 



other utilized a frequency comb based on a fiber laser 
located at California State University - East Bay. 



A. The Ti:Sapphire-Comb Experiment 

The frequency comb based on the Ti:Sapphire mode- 
locked laser is similar to the "Standard Ti:Sapphire" res- 
onators discussed in Ref. [13] ■ The mirrors of the cav- 
ity consist of two mirrors with a radius of curvature of 
30 cm and two plane mirrors. The group- velocity dis- 
persion (GVD) of the curved mirrors is compensated 
with a net negative group- velocity dispersion of —100 fs 2 
over a bandwidth of 620-1050 nm. One of the plane 
mirrors also has negative group velocity dispersion of 
—40 fs 2 . The other plane mirror is 99% reflective and 
serves as the output coupler. The Ti:Sapphire crystal 
has an absorbance of a — 5.65 at 514 nm and an optical 
path length of 1.5 mm. The repetition rate of the laser 
is w 925 MHz. The TkSapphire resonator is pumped 
with 5.5 W of light from a diode-pumped frequency- 
doubled neodymium vanadate laser at 532 nm. The spec- 
tral bandwidth of the output of the TkSapphire laser is 
ss 30 nm, centered at w 780 nm with an output power of 

550 mW. This light is coupled into ~ 30 cm of a non- 
linear photonic crystal fiber with a zero GVD wavelength 
at 790 nm. The light is spectrally broadened by the fiber 
to span rj 500-1100 nm. Approximately 200 mW of light 
is coupled through the photonic crystal fiber. The output 
of the fiber is split according to its spectral region by use 
of a short-pass interference filter. The spectral regions 
from 500-650 nm and from 1000-1100 nm are used for 
the stabilization of the carrier envelop offset frequency, 
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fo, and the repetition rate, f r , of the optical frequency 
comb, while the remaining spectral region can used for 
spectroscopy. 

The carrier envelope offset frequency is detected via 
the self-referencing technique described in Ref . [H] . The 
infrared light near 1060 nm is doubled using a peri- 
odically poled lithium niobate crystal with a length of 
1 mm that is thermally stabilized to 381 K. The dou- 
bled infrared light and the green light that is directly 
produced in the fiber are filtered by use of an inter- 
ference filter centered at 530 nm with a pass band of 
30 nm and detected with an amplified photodiode with 
a bandwidth of 2 GHz. The interference of the doubled 
infrared light with the directly generated green light re- 
sults in a radio-frequency signal at the carrier-envelope 
offset frequency. The signal-to-noise ratio of fo is typ- 
ically 40 dB in a 300 kHz resolution bandwidth. The 
carrier-envelop offset frequency is phase-locked to a signal 
generator using standard phase-locking techniques (see 
e.g. Ref. [H]). The line width of the offset frequency is 
typically < 500 kHz. The repetition rate is detected with 
the same detector that is used to detect fo. The output 
of the photodiode is split and filtered to select the rep- 
etition rate. The repetition rate is phase stabilized to 
a second signal generator using an rf mixer. The sig- 
nal generators are referenced to a rubidium atomic clock 
that is steered to the global position system (GPS) re- 
sulting in a fractional frequency uncertainty of ~ 10~ n 
in 1 second. 

The 650-1000 nm light generated from the nonlinear 
photonic crystal fiber is passed through an interference 
filter to select the wavelengths of interest. For excitation 
through the 5P3/2 intermediate state the light is passed 
through an interference filter centered at 780 nm with 
a 20 nm pass band. For excitation through the 5Pi/ 2 
state the light is passed through an interference filter 
that transmitted light from ps 785-860 nm. This spectral 
range allows for excitation of the 5S;iy 2 — > SP1/2 5D 3 / 2 
and the 5Si/2 — > 5P 3 / 2 — > 5D 3 / 25 / 2 transitions. The 
light is focused with a 10-cm lens onto a rubidium vapor 
cell at s» 328 K. The vapor cell is a cylinder with a 25 mm 
diameter and a 75 mm length. The light is propagated 
through the curved wall of the vapor cell. A second lens 
with a focal length of 10 cm is used to re-collimated the 
laser beam. The light is retroreflected off a mirror and 
sent back through the vapor cell. A small fraction of the 
light is picked off for power monitoring by use of a 8% 
reflective pellicle beam splitter. 

Atoms excited to the 5D 5 / 2 state decay via cascade 
through the 6Pi/2 l 3/2 states, emitting light at 420 nm. 
This fluorescence is detected with a photomultiplier tube 
with an active area of 0.50cm 2 . An interference filter cen- 
tered at 420 nm with a pass band of 20 nm and a short- 
pass filter at 650 nm are placed in front of the photomulti- 
plier tube to reduce background and scattered light. The 
photocurrent from the photomultiplier tube is sent to a 
transimpedance amplifier with a gain of 10 V/A and 
a bandwidth of 20 kHz. The output of the amplifier is 



digitized with a analog-to-digital converter and recorded 
on a computer. The repetition rate of the optical fre- 
quency comb is stepped in 1 Hz intervals over a region of 
« 1.5 kHz and the fluorescence is recorded for ps 1 second 
at each frequency step. 



B. The Fiber-Comb Experiment 

The frequency comb based on the erbium-doped fiber 
laser is a commercial system from Menlo Systems GmbH 
(FC1500-250-WG Optical Frequency Synthesizer) con- 
sisting of a femtosecond laser system (M-Comb) with 
repetition rate f r ps 250 MHz whose output is con- 
nected via fiber-optic cables to (1) a fast PIN photodi- 
ode that detects the fourth harmonic of f r (for stabiliza- 
tion/control of / r ), (2) an erbium-doped fiber amplifier 
(M-Phase EDFA) whose output is sent to an / — 2/ inter- 
ferometer (XPS-WG 1500) used to measure the carrier- 
envelope offset frequency f (for stabilization/control of 
fo), and (3) a second erbium-doped fiber amplifier whose 
output is directed through a second-harmonic generating 
(SHG) crystal. The output power after the SHG crys- 
tal is typically ps 230 mW centered about 780 nm with 
a FWHM bandwidth of ps 15 nm. The laser light is 
then focused into a photonic crystal fiber (PCF) which 
spectrally broadens the output to span roughly from 600- 
900 nm. The power after the PCF is typically ps 85 mW. 

The fourth harmonic of f r detected by the fast PIN 
photodiode is mixed with a 980 MHz signal generated 
by a dielectric resonator oscillator (DRO) to produce a 
20 MHz intermediate frequency, which is counted and 
further mixed with the ps 20 MHz output of a tunable low 
frequency direct-digital synthesizer (DDS), which enables 
adjustment of f r around its nominal value of 250 MHz. 
The DC signal resulting from the second stage of mixing 
is used in a feedback loop to lock f r by controlling an 
intra-cavity piezo in the M-Comb laser. All synthesizers 
and counters are referenced to a 10 MHz signal from a 
GPS-disciplined, ultrastable quartz oscillator (TimeTech 
Reference Generator) with relative stability better than 
5 x 10~ 12 in 1 second. 

The offset frequency fo, which is tuned to 20 MHz, 
is counted and fed into a lock-in detector. The latter is 
referenced to a 20 MHz signal directly generated from the 
clock signal by frequency doubling. The lock-in detector 
output signal is used in a feedback loop to lock fo by 
controlling the M-Comb laser pump power. 

The laser light spanning 600-900 nm generated by the 
PCF is then split into two beams with a 50/50 nonpolar- 
izing beamsplitting cube. The beams pass through 10- 
nm bandpass interference filters centered near the wave- 
lengths of the transitions of interest. One of the beams 
(corresponding to the second stage of the two-photon 
transition of interest) passes through a chopper wheel 
which modulates the light at ss 250 Hz. Mirrors direct 
the two beams so that they are counterpropagating and 
are focused by antireflection-coated, 10-cm focal length 
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bi-convex lenses at the center of a cylindrical Pyrex cell 
(25 mm in diameter and 75 mm in length) containing 
Rb vapor heated to f=s 323 K. The beams enter per- 
pendicular to the flat circular cell windows and counter- 
propagate along the axis of the cell. The focal point 
of the beams is at the center of the cell. Fluorescence 
from the cascade decay of the upper 5D 5 / 2 and 5D 3 / 2 
states at 420 nm is monitored by an amplified photomul- 
tiplier tube (ThorLABs PMM01, active area diameter 
= 22 mm) fitted with an interference filter centered at 
420 nm with a 10 nm pass band to reduce background 
light. The PMT is located to the side of the cell di- 
rectly viewing the center of the cell where the counter- 
propagating beams are focused. Approximate alignment 
of the counter-propagating beams is achieved by directing 
the two focused beams through 100 fim pinholes without 
the cell in place. The alignment is subsequently improved 
by adjusting the position of the lenses to maximize the 
detected fluorescence signal when f r is tuned to a res- 
onant value for the two-photon transition. The output 
of the amplified PMT is sent to the input of a lock-in 
amplifier (Signal Recovery Model SR7265) referenced to 
the chopper wheel. The demodulated output of the lock- 
in amplifier (time constant = 2 s) is recorded as f r is 
stepped in 0.25 Hz increments every 4 s. 



IV. RESULTS AND DISCUSSION 
A. The Ti:Sapphire Comb Data and Calculations 

The experimental and calculated spectra for the 
5Si/ 2 -> 5Pi/ 2 -> 5D 3/2 and 5Si/ 2 -> 5P 3 / 2 -> 5D 5/2 
transitions from the Ti:Sapphire-comb experiment are 
shown in Figs. [4] and [5] The data in Fig. [4] were taken 
with an interference filter that had a stronger transmis- 
sion at 794 nm (resonant wavelength for the 5Si/ 2 —> 
5Pi/2 transition) than at 780 nm (resonant wavelength 
for the 5S — > 5P 3 / 2 transition). This results in a spec- 
trum largely consisting of the 5Si/ 2 — > 5Pi/ 2 — » 5D 3 / 2 
transitions, with much smaller contributions from the 
5Si/2 — > 5P 3 / 2 — » 5D 3 / 2 5 / 2 . transitions. The data in 
Fig. [5] were taken with an interference filter that did not 
transmit light at 794 nm so that the spectrum consists 
entirely of the 5Si/ 2 — > 5P 3 / 2 — ► 5D 3 / 2 ,5/2 transitions. 
There is excellent agreement in the peak positions for 
both sets of data. The 5Si/ 2 — > 5P 3 / 2 — > 5D 3 / 2 transi- 
tions are much smaller than the 5Si/ 2 — » 5P 3 / 2 — » 5D 5 / 2 
transitions due to angular momentum considerations. 

The amplitude agreement is reasonable, given the as- 
sumptions of the model. We attribute the residual dis- 
agreement between the calculated and experimental am- 
plitudes to issues related to the power depletion of the 
light exciting the first stage of the transition. The cal- 
culation includes this effect but assumes that the coun- 
terpropagating beams have common focal points which 
localizes the position at which the atoms are excited. If 
the focal positions of the two laser beams do not overlap 



perfectly, the location where the atoms are excited can 
vary depending on the relative powers of the two beams. 
This can combine with position-dependent fluorescence 
detection efficiency and affect the relative amplitudes of 
resonances from different velocity classes. We note that 
an effective optical depth of 1.0 cm was used in the cal- 
culation, which agrees well with geometrical size of the 
cell. 

The calculations were done with a transit line width 
of tt 1 = 4 MHz. This value is somewhat larger than the 
estimated line width if If- = 500 kHz, based on a laser 
beam diameter of ~ 80 /im. We attribute this additional 
broadening to imperfect overlap in the focal points of the 
counter-propagating beams. If the foci of the two beams 
do not overlap completely then there will be atoms for 
which the wave vectors of the two beams are not an- 
tiparallel. These atoms will have an additional broaden- 
ing due to the imperfect cancelation of the Doppler shift. 
We explored the effect of wave vectors that were slightly 
misaligned by including a misalignment term in the cal- 
culation. The calculation reproduced broadening at the 
same level we observe in the experiment for values that 
are consistent with our limits of how well the foci of the 
two beams overlap. The calculations were done assum- 
ing linearly polarized light. For the data, the polarization 
of the the light was not controlled. However, data were 
taken with linearly polarized light and no differences in 
the spectra were observed. 

Trace (c) of Fig. [4] shows the calculated spectrum for 
5Sj./2 -> 5Pi/2 -> 5D 3 / 2 transition for 85 Rb. The spec- 
trum consists of three distinct groups. The first group 
on the left is due to excitation from the F — 3 — > F" 
transitions, while the second grouping corresponds to the 
F = 2 — > F" transitions. The last group at the high 
end of the spectrum corresponds again to the excitation 
of the F = 3 — > F" transitions, but with a new total 
mode number equal to one fewer than the mode number 
corresponding to the first group. The frequency spacing 
between the first and third group corresponds the the 
repetition interval predicted by Eq. Af r w 1.1 kHz. 
The relative amplitudes of the groups, as well as the am- 
plitudes of the peaks within each group, vary significantly 
due to the different velocity classes that are excited by 
the different modes of the frequency comb. Such features 
are also evident in the 5Si/ 2 — > 5P 3 / 2 — > 5D 5 / 2 transi- 
tions (Fig. [5]). 

In both spectra the hyperfine states of the ground and 
excited states are clearly resolved while the intermediate 
hyperfine structure is not. This is in contrast to the 
data in Ref. [6] where the transitions through different 
intermediate states gave distinct peaks. This difference is 
a result of the relatively close energies of the two resonant 
photons. 
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FIG. 4: Color Online: Trace (a) is the experimental spectrum 
taken with the TkSapphire comb. Trace (b) shows the cal- 
culated fluorescence spectra for the 5Si/2 — > 5Pi/2 — > 5D 3 / 2 
and 5Si/2 — > 5P 3 /2 —¥ 5D 3 / 2 ,5/2 transitions. Traces (c) and 



(d) show the calculation for the 5Si 
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FIG. 5: Color Online: Trace (a) is the experimental spectrum 
taken with the Ti:Sapphire comb. Trace (b) shows the calcu- 
lated fluorescence spectra for the 5Si/2 ~~ ► 5P3/2 — > 5D 5/ / 2 - 
Traces (c) and (d) show the calculation for the 5Si/2 — > 
5P3/2 5D 5 /2 transitions for 85 Rb and 87 Rb, respectively. 
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transitions for 85 Rb and 87 Rb, respectively. Trace (e) shows Trace (e) shows the calculation for the 5Si 
the calculation for the 5Si 
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transition 5D 3 / 2 transition for both isotopes. This transition is nine 



for both isotopes. The amplitude of the transitions through 
the 5P3/2 state relative to the those through the 5Pi/ 2 state 
were adjusted to reflect the differing light intensity transmit- 
ted through the filter. 



-"3/2 

times weaker due to the relative amplitudes of the reduced 
matrix elements and is also weighted by the branching ratio 
to the 6P states (38% for the 5D 3 / 2 state and 35% for the 
5D5/2 state based on the calculations of Ref. [16) V 



B. The Fiber Comb Data and Calculations 

The experimental and calculated spectra for the 
5Si/ 2 -> 5Pi/ 2 -> 5D 3/2 and 5Si/ 2 -> 5P 3/2 -> 5D 3/2i5/2 
transitions from the fiber-comb experiment are shown in 
Figs. [Bland [7j respectively. For the spectra shown in Figs. 
[6]andj7j interference filters with 10- nm bandwidths were 
employed. These interference filters provided sufficient 
wavelength selectivity so that in Fig. [6] only excitation 
through the 5Sx/2 ~> 5Pi/ 2 pathway is observed and in 
Fig. [7| only excitation through the 5Si/2 — > 5P 3 / 2 is ob- 
served. 

As with the Ti:Sapphire-comb experiment, the rel- 
ative amplitudes of the different resonance peaks dif- 
fer slightly between the calculated and experimental 
spectra. The position of the peaks agree well for the 
5Si/2 — > 5P 3 / 2 — > 5D 5 / 2 spectrum. However, for the 
5Si/2 — > 5Pi/ 2 — > 5D 3 / 2 spectrum there is some dis- 
crepancy at the 0.5 Hz level. We attribute this to the 
increased overlap and complexity of the spectrum. Shifts 
in the apparent positions of the transition resonances will 
occur as a result of the discrepancies in relative ampli- 
tudes between the calculated and experimental spectra. 

Power depletion of the light corresponding to excita- 
tion of the first stage of the transition was not found 



to improve the agreement between the data and calcu- 
lation and was not included. This may be a result of a 
lower atomic density for the data taken with the fiber 
comb compared to that for the ThSapphire comb exper- 
iment. The calculations for the 5Si/ 2 — > 5Pi/ 2 — > 5D 3 / 2 
transitions were done with a transit line width of 7^ = 
2.7 MHz. Those for the 5S 1/2 -> 5P 1/2 -> 5D 3/2 transi- 
tions were done with a transit line width of 5^ = 4 MHz. 
These values are in reasonable agreement with the es- 
timated based on a laser beam diameter of 15 /im, al- 
though additional broadening due to imperfect overlap 
of the foci of the counter propagating laser beams may 
have contributed additional broadening. The calcula- 
tions were done assuming linearly polarized light. For 
the data shown, the polarization of the the light was not 
controlled. However, data were taken with linearly po- 
larized light for the transitions excited through the 5P 3 / 2 
state and these data did not display any differences from 
the data where the polarization was not controlled. 

For a comb with a repetition rate of 250 MHz, f r must 
be scanned by Af r rs 80 Hz in order for a given transition 
to recur in the spectrum. As can be seen from Eq. (|7|, 
this is a factor of the ratio of the repetition rates squared 
smaller than the repetition interval for the ThSapphire 
comb. The widths of the resonance peaks decrease by 
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a single power of the ratio of the repetition rates. As 
a result, there is more overlap of the peaks arising from 
different transitions. This is clearly observable in the 
comparison of 5S]/2 — > 5P 3 / 2 — > 5D 5 / 2 transitions for 
the different combs, Figs. [5] and [7} However, in the case 
of the 5Si/2 — >• 5P3/2 — > 5D5/2 transition the individual 
peaks are still resolvable in most cases. Trace (f) shows 
the calculated spectrum for a single hyperfine transition, 
in this case the F = 2 -t F" = 3 for 85 Rb. As expected, 
there is only one resonant peak in the spectrum for this 
transition, a result of the near degeneracy of the two 
photons involved in the transition. 

As described above, the situation is quite different for 
the 5Si/2 — > 5Pi/2 — > 5D 3 /2 transitions. Trace (e) of 
Fig. [6] shows the calculated spectrum for a single hyper- 
fine transition of 85 Rb (F = 3 ->• F" = 2). There are 
multiple distinct velocity classes that contribute to the 
fluorescence signal and the resonance for each velocity 
class occurs at a different repetition rate. The result is 
multiple distinct peaks in the fluorescence signal. This 
results in a significantly more complicated fluorescence 
spectrum where there are numerous peaks overlapping 
peaks. 
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FIG. 6: Color Online: Trace (a) is the experimental spectrum 
taken with the fiber comb. Trace (b) shows the calculated 
fluorescence spectra for the 5Sx/a ~~ ► 5Pi/2 — > 5D 3 / 2 transi- 
tions. Traces (c) and (d) show the calculation for the 5Si/2 — >• 
5Pi/2 ~ > 5D3/2 transitions for 85 Rb and S7 Rb, respectively. 
Trace (e) shows the calculation for the F = 3 — > F" = 2 tran- 
sition for 85 Rb. Several distinct resonances occur for different 
mode combinations, each of which have the same total mode 
number. 
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FIG. 7: Color Online: Trace (a) is the experimental spec- 
trum taken with the fiber comb. Trace (b) shows the calcu- 
lated fluorescence spectra for the 5Si/2 — > 5P3/2 — ¥ 5D 3/ / 2 ,5/2 
transitions. Traces (c) and (d) show the calculation for the 
5Si/ 2 — > 5P 3 /2 — > 5D 5 /2 transitions for 85 Rb and 87 Rb, re- 
spectively. Trace (e) shows the calculation for the 5Si/2 — > 
5P3/2 5D3/2 transition for both isotopes. This transition 
is nine times weaker due to the relative amplitudes of the re- 
duced matrix elements and is also weighted by the branching 
ratio to the 6P states (38% for the 5D 3 / 2 state and 35% for 
the 5D5/2 state based on the calculations of Ref. [16]). Trace 
(f) shows the calculation for the F — 2 — > F" = 3 transition 
for 85 Rb. 



V. CONCLUSIONS 

We have presented a comprehensive investigation of 
velocity-selective two-photon direct frequency comb spec- 
troscopy in atomic Rb. We have experimentally and the- 
oretically demonstrated the effect of the repetition rate 
of a frequency comb on the two-photon excitation rate. 
The energy level structure of Rb also allowed us to ex- 
plore the effect of the energy of the intermediate state 
on the two-photon excitation. The energy of the inter- 
mediate state is particularly important in the case where 
the repetition rate of the frequency comb is less than the 
Doppler width of the resonance corresponding to the first 
stage of the transition. 

This investigation demonstrates the benefits and chal- 
lenges of direct-frequency-comb spectroscopy. While 
there is significant advantage in the wavelength versatil- 
ity of the frequency comb, the details of resulting spectra 
can display complicated features resulting from the pres- 
ence of the numerous frequencies. This effect is partic- 
ularly pronounced in the case where the repetition rate 
of the comb is less than the Doppler distribution of the 
atoms. The presence of multiple resonances correspond- 
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ing to excitation of a two-photon resonance from comb 
light with different mode numbers results in a dense spec- 
trum. This effect will likely complicate any effort to make 
quantitative measurements of transition frequencies for 
atoms with a multiple hypcrfine transitions. 
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